Identification of strong and weak interacting two level systems in KBr:CN 
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Tunneling two level systems (TLSs) are believed to be the source of phenomena such as the uni- 
versal low temperature properties in disordered and amorphous solids, and 1/f noise. The existence 
of these phenomena in a large variety of dissimilar physical systems testifies for the universal nature 
of the TLSs, which however, is not yet known. Following a recent suggestion that attributes the low 
temperature TLSs to inversion pairs [M. Schechter and P. C. E. Stamp, arXiv:0910.1283 ] we cal- 
culate explicitly the TLS-phonon coupling of inversion symmetric and asymmetric TLSs in a given 
disordered crystal. Our work (a) estimates parameters that support the theory in M. Schechter and 
P. C. E. Stamp, arXiv:0910.1283 in its general form, and (b) positively identifies, for the first time, 
the relevant TLSs in a given system. 



Introduction — Amorphous solids and many disor- 
dered lattices show peculiar universal characteristics at 
low temperatures [l|-0]- Below Tjj « 3K systems which 
are otherwise very different have specific heat C v oc T a , 
with awl, thermal conductivity k oc T' 3 with (3 « 2, and 
internal friction Q m 2ttI/X w 10 3 , independent of T, A, 
and with only a small variance between materials. Here I 
is the phonon mean free path and A is the phonon wave- 
length. In an effort to explain this remarkable universal- 
ity, Anderson Halperin and VarmaQ and Philips^ sug- 
gested a phcnomenological theory, where the existence 
of tunneling two level systems (TLSs) in these materials 
was postulated, and an ansatz for their density of states 
was given. This "standard tunneling model" (STM) has 
been very successful in explaining the above mentioned 
phenomena. Still, the identity of the tunneling TLSs has 
remained unknown. Furthermore, the smallness and uni- 
versality of the phonon attenuation, and the energy scale 
dictating Tjj are not accounted for by the STM. 

Two level systems are also believed to be the cause of 
1/f noise. Recently, it has been shown that 1/f noise 
is the main source for decoherence of superconducting 
qubits, and a major obstacle in their ability to perform 
quantum computation^. Also in these systems the na- 
ture of the TLSs is not known, yet assuming their exis- 
tence and applying the STM has resulted in an explana- 
tion of the low frequency 1/f noise and high frequency 
linear in / noise on the same footing Q. 

Extensive experimental investigations have revealed 
that the condition to observe universality is the pres- 
ence of tunneling states and strong lattice strain [H, Q, 
and that the phenomena in amorphous solids and dis- 
ordered crystals are equivalent Disordered crystals 
are advantageous for both experimental and theoretical 
investigation (IT) . Experimentally, they allow control of 
the nature and relative concentration of host material 
and impurities, and therefore a detailed study of differ- 
ent universal properties and their origin. The existence of 
lattice structure and the apparent candidates for tunnel- 
ing states allows a favorable starting point for theoretical 
treatment as well. 



Indeed, it was argued [12| that, at least in disordered 



crystals, tunneling states can be categorized into two 
types of TLSs, denoted r and S. The states of a r-TLS 
are related to each other by inversion. Consequently, the 
interaction of a r-TLS with the phonon field 7w is small, 
as it results only from disorder induced local deviations 
from inversion symmetry. The S'-TLSs are asymmetric 
with respect to local inversion, with a stron g in teraction 
with the phonon field j s . It was then shown [l2| that the 
S'-TLSs are gap ped below Tjj by the r-TLSs through an 
Efros Shklovskii[l3| type mechanism, and that below Tjj 
the r-TLSs arc effectively noninteracting, and dictate the 
phonon attenuation. Thus, at T < Tjj the r-TLSs fulfil 
the assumptions of the STM. The small parameter of the 
theory is g = 7 w /7s ~ E<p/E c ~ (1 — 3) x 10~ 2 , where 
, Ec are the typical elastic and Coulomb energies in 
the system. This small parameter gives the universality 
and smallness of the phonon attenuation. Defining Tq 
as the ordering temperature of the S'-TLSs, the emerg- 
ing DOS of the S-TLSs at an energy Tjj fts gT G , dic- 
tates Tjj as the energy scale below which universality is 
observed fl2|. 

In this Letter we use DFT and ab-initio calculations to 
calculate the interaction of TLSs of types r and S with 
the phonon field in the system KBr( 1 _ x )(CN) :r (KBr:CN, 
Fig. EE). Wc find that 7w w O.lcV, and 7 S rj 3eV. Our 
estimation of 7 W compares well with the experimentally 
measured value for the relevant TLSs at low energies, 
of 7 0.1 2cV for im pur ity concentration x= 0.25 and 
7 0.2cV for x= 0.5[1J, Qjl]. Our results also support 
the central arguments of the theory in Ref.[l2[ in (i) the 
categorization of the TLSs according to their symmetry 
under inversion, (ii) the ratio of the strengths of their in- 
teractions with the phonon field, constituting the small 
parameter of the theory, and (iii) the identification of the 
symmetric TLSs as the relevant TLSs dictating the low 
temperature universal properties in disordered solids. In 
addition, we re-enforce the prediction made in Ref.[l2l| 
for the existence, at higher energies, of a second type of 
TLSs (of S type), with a much stronger coupling to the 
phonon field. Note, that although we focus here on the 
simplest single impurity excitations, our analysis does 
not exclude the possibility of symmetric and asymmetric 
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multi- impurity excitations [16M18I ] . Such excitations are 
expected to be significant especially for systems where 
single impurity excitations do not produce symmetric 
TLSspSlIj. 

For the specific KBr:CN system, early works have sug- 
gested, based on theories quite different from that of 
Ref . 11211. that CN flips comprise the relevant low energy 
TLSs|2ll. |22| . However, for long, advance in this direc- 
tion was hindered because of experiments showing that 
the substitution of the symmetric N 2 molecules for the 
asymmetric CO molecules in N ? /A r/CO does not change 
its universal characteristics [l9l [20| . Our results here, in 
conjunction with the theory in Ref. (l2"| . positively iden- 
tify the 180° CN flips as the relevant TLS excitations dic- 
tating the low temperature characteristics in the KBr:CN 
system. Reconciliation of our results with the experi- 
ment in Ref.pjl, [2(| stems from the fact that pairs of 
N2 molecules do produce symmetric TLSs in the ArN 2 
system [2^|. 
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FIG. 1: Three 5x5 fragments of a KBr:CN lattice. The up 
state of the central impurity in fragment (b) is related by a 
r excitation to a down state in fragment (a), and by an S 
excitation to a left state in fragment (c). 

Calculation — KBr:CN is perhaps the most studied 
disordered lattice showing universal characteristics. The 
CN - impurities have been found to orient either in the 
direction of the in-space diagonals, preferred for very low 
CN - concentrations 
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tions at high temperatures [25| . or in the direction of the 
axes, preferred for intermediate CN - concentrations at 
low temperatures [1^. The six (eight) possible states of 
each impurity can be categorized into three (four) inver- 
sion pairs, each having two states related to each other by 
an 180° flip. Such flips constitute r excitations, whereas 
rotations between different axis (diagonals) correspond 
to S excitations (l2l [26| . 

The interaction of such a system with the lattice can 
be described by the Hamiltonianjl^, [26| 

H ** = EE^ + ^ S i + -tifrf] (1) 

3 a,0 

where 77 is an orientation-independent volume factor and 
Uapi^j) denotes the phonon field at point rj. Whereas 
the central purpose of this Letter is the calculation of 
7 W and j s , we also calculate the parameter 77 for both 
CN - and CI - impurities. This parameter determines 
the strain, and thus the effective random field in the 
system [2^, [27j. Usually 77 < 7 S . In KBr:CN 77 is signifi- 
cantly subdominant, as the Br - and CN - ions have sim- 
ilar volumes 15 1. In KBr:Cl this term is responsible for 



the strains allowing for the existence of universal proper- 
ties upon minimal CN - dilution^, 0- This random field 
term is also central to the smearing of the glass transition 
and the peculiar disordering of dilute glasses 27j . 

Following the above definition of 77, 7 S and 7 W , we de- 
vise a series of numerical calculations to estimate them. 
In sum, we choose a number of lattice fragments and 
use DFT/ab initio methods to calculate the energy dif- 
ference between the effects of phonon-likc perturbations 
of the system with a central CN - impurity in different 
states. For simplicity and without loss of generality, we 
restrict the excitations of the CN - impurity under study 
to two dimensions, thus the possible states are up, down, 
left and right. To those possible orientations, we apply 
vertical or horizontal phonons. Usually the symmetry is 
low enough to allow for several independent estimations 
for each parameter. 

Our determination of 7 W and 7 S is performed as follows: 



7w,s = i- 1(^(6)--^) 



iE 3 ph {b) 
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where {i, j} are {up, down} or {left, right} for 7 W and 
{up, left} or {down, right} for 7 S , ph can stand for 
vertical or horizontal phonons, E l is the energy of an 
impurity i surrounded by a lattice fragment in its equi- 
librium geometry and E l h is the energy of the same im- 
purity after a lattice contraction by a fraction b along 
a given crystallographic coordinate, mimicking the effect 
of a longitudinal phonon. For 77, the same procedure is 
applied where {i,j} means presence or absence of impu- 
rity, and for CN - impurities all possible orientations are 
averaged. 
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(a) (b) (c) 

FIG. 2: Atomic coordinates of a lattice fragment in three key 
steps of a sample calculation, (a): all positions are optimized, 
without a central TLS, but in presence of an impurity (2,0,u). 
(b): a TLS (0,0, u) is substituted in the place of the central 
Br - , its position is optimized, yielding E % . (c): the lattice is 
contracted by 2% vertically, yielding El ertical (— 2%). 

Fig. [2] illustrates the detailed procedure of a 7 calcu- 
lation. We describe the position and orientation of the 
impurity as (2,0,u), where the numbers refer to the coor- 
dinates and the letter to the orientation, in this case the 
CN - is lying on the abscise, at 2 interatomic spacings 
to the right, and the nitrogen is pointing up. We then 
proceed with the following steps: (i) for that particular 
system, the atomic positions are found which minimize 
the energy in the absence of the TLS under investiga- 
tion, i.e., when the central position is occupied by a Br - ; 
that is our definition of the equilibrium geometry of the 
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lattice, (ii) a CN~ (facing "up", in this case) is substi- 
tuted for the central Br~, and, freezing the lattice, only 
the position of the CN~ is optimized, to obtain E up . 
(iii) the atomic positions are contracted along the verti- 
cal axis to obtain , (6) . The calculation is done for 
b E {-5%, -2%, -1%, -0.25%, 0.25%, 1%, 2%, 5%}, and 
the limit of small b is taken (see Fig. [3]). The repetition of 
this procedure for the vertical phonon and the " down" , 
"left" and "right" orientations of the central CN~ allows 
for four non-independent estimations of both 7 W and 7 S 
(in 3D wc also have the "front" and "back" orientations). 
In the case of r), the weighted average for all orientations 
yields one unique estimation. Note that in step (ii) we 
do not relax the whole lattice in the presence of the TLS. 
Our procedure is in line with both processes of lattice 
relaxation and phonon scattering by TLSs, which result 
from the out of equilibrium first order interaction of the 
TLS with the lattice. 



second-order perturbations. Last, we include a limited 
study of larger samples with a higher number of impuri- 
ties, using HFS calculations with the minimal basis STO- 
3G on a 7x7 fragment with up to 4 extra impurities. 




FIG. 3: Energy response E (in meV) to horizontal phonons of 
amplitude U. 7 W (main graph) and 7 S (upper inset) as linear 
fits on selected fragments on table [T] open circles: (b); filled 
circles: (c); squares: (d). Lower inset: second-order fit for n 
(see text); open circles: Ear — E up ; filled circles: Ebt — Ei e f t ; 
squares: weighted average for all orientations. 



We use the standard package Gaussian03 [28[ to per- 
form quantum chemistry calculations on lattice frag- 
ments of different sizes and shapes and at different lev- 
els of sophistication. As we are dealing with a local 
phenomenon, and for cost reasons, most of the calcula- 
tions are performed on small zero-dimensional squares or 
cubes, either 3x3, 3x3x3 or 5x5. The TLS under evalua- 
tion is always in the center, so that any deviation from a 
ccntrosymmetric situation felt by the TLS is due to the 
extra impurities and not to border effects. We mainly 
use the hybrid DFT/ ab initio method B3LYP with small 
orbital sets, either 3-21G or 6-31G. The influence of a 
better description of the anions is tested by repeating 
some calculations with the more flexible basis sets, up 
to 6-311+G*. Additionally, we check the relevance of 
dinamical correlation by comparing plain Hartrce-Fock 
with the Moeller-Plesset perturbation theory to the sec- 
ond order (MP2), which includes double excitations as 



fragment 


impurity 


method 


7s 


7™ 


3x3 


none 


HF/6-311+G* 


3.02 


(sym) 


3x3 


none 


MP2/3-21G 


1.91 


(sym) 


3x3 


none 


MP2/6-311+G* 


3.55 


(sym) 


3x3 


none 


B3LYP/3-21G 


3.28 


(sym) 


3x3 (tr.ph.) 


none 


B3LYP/3-21G 


2.95 


(sym) 


3x3x3 


none 


B3LYP/3-21G 


4.20 


(sym) 


3x3 (tilt) 


(1,1, dr) 


B3LYP/6-31G 


3.40 


0.08 


3x3<-> 


(l,-l,dl)(l,l,ul) 


B3LYP/3-21G 


2.86 


0.03 


3x3<"> 


(l,-l,dr)(l,l,ul) 


MP2/6-311G 


1.80 


0.003 


3x3x3 (c) 


(0,-l,l,r)(l,-l,0,f) 


HF/6-31+G 


2.15 


0.04 


3x3x3 (tilt) 


(l.l.dl) 


B3LYP/3-21G 


n.a. 


0.14 


5x5< d) 


(0,-2,r) 


B3LYP/6-31G 


4.70 


0.04 


5x5 




B3LYP/3-21G 


2.10 


0.04 


5x5 


(2,0,d) 


B3LYP/3-21G 


2.40 


0.11 


7x7<=> 


(see caption) 


HFS/STO-3G 


10 - 20 


0.1 - 0.2 



TABLE I: Some estimations, in eV, for 7 S and 7 W . Absolute 
values are given, since the TLS orientation is arbitrary. Ex- 
cept for fragments (b-d), all possible TLSs were calculated, 
as shown for fragment (a) in table [TTJ and only the highest 
values are displayed here. In the absence of extra impuri- 
ties, 7 W = for symmetry reasons. For fragments (b-d), only 
one r-TLS and one S-TLS were chosen, and linearity of the 
energy response was tested as shown in Fig. [3] For these frag- 
ments the values for 7 W should be considered as lower bounds, 
tr.ph. stands for "transverse phonon". Fragment (e) sums up 
three calculations with a minimal basis set, with up to four 
impurities at positions (0,-2,r)(-2,2,r)(2,0,u),(2,2,d); these re- 
sults should not be taken on equal footing with the rest of the 
table. 

Results — Fig. [3] illustrates some tests of the range 
of linearity. One can see that the results are essentially 
the same for different fragments and levels of calcula- 
tions: the first order approximation is very accurate at 
least for phonon amplitudes of one or two percent of the 
interatomic spacing. On a 3x3 fragment with a central 
CN" impurity, at MP2/6-31+G level, Fig. |3] shows the 
estimation of r\ = O.QeV. The same conditions yield a 
comparable rj = 0.9eV" for a Cl~ impurity. In all cases 
a noticeable second-order correction of the order of 5eV 
can be fitted. 

The central result of this Letter, reported in table HI is 
the calculation of 7 S ~ 3eV and < j w < 0.15eV. The 
finite size of our samples, the quality of our calculation 
methods, and differences between planar and cubic sam- 
ples, all lead to some variance in the parameters. Yet, 
the strength of our results lies in the fact that our esti- 
mations are fairly consistent in their order of magnitude 
for very different lattice fragments and a variety of levels 
of calculation. This is true for additional calculations, 
e.g. for a non-central CN~ impurity, not reported here. 
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7= 






7w 


A ul 


2.84 




Ah d 


0.0036 


Af 


2.82 




a5; 


0.0070 


K l 


2.86 






0.0001 


Af 


2.83 




A (r 


0.0273 



TABLE II: Values in eV, summary of results for fragment 
(a) on Table ffl Shorthand Af = (EZ rlzontal - E^) - 
(EfcTzontal - E down ) has been used for clarity. 

The calculations using a minimal basis set serve to dis- 
card a correlation between j w or 7 S and the number of 
impurities. 

One should point out that within a given sample and a 
given level of calculation, the variance in the values of 7 S 
between different orientations is a result of the small elas- 
tic deviations from symmetry, and are therefore a factor 
of g smaller than its typical value. With regard to j w , 
its values are dictated by the aforementioned deviation 
from local inversion symmetry. Thus, the distribution of 
all possible estimations of its values is peaked at zero, 
with a variance which equals the typical value. This is 
displayed in table HH for all TLSs in fragment (a). It 
can also be seen in the main panel of Fig. [31 where the 
particular r-TLS chosen for fragment (b) experiences a 
very symmetric environment. Note that for each particu- 
lar combination of fragment and calculation method, the 
highest 7 W values obtained among the four TLS-phonon 
combinations are denoted as our estimate values for 7 W 
in Table HI These values are the most relevant for our 
purposes, as they are expected to be the best predictors 
for a real system with many impurities. 

A noteworthy complication is presented by the in-plane 
diagonal orientations of the TLS, because of the small 
energy difference with the in-space diagonal states in the 
real system. Depending on the fragment, impurities and 



calculation method, the relative order of stability changes 
and the energy of one or more of the " axial" orientations 
rises above the most stable "diagonal" orientation. Some 
examples of these orientations, denoted as {dl, ul, dr, 
ur}, are shown in table [I] As illustrated in Figs. [1] and [2] 
in the perimeter of the fragment, where the CN _ suffer 
from intense border effects, we even find intermediate ori- 
entations. In the cases where the central TLS is affected 
by this problem -noted in Table U as "tilt"- the extraction 
of the parameters can be technically more difficult, but 
there is no fundamental physical difference between axial 
and diagonal orientations as in both cases there is a clear 
distinction between S-TLSs and r-TLSs. 

Summary — Our numerical calculations confirm qual- 
itatively and quantitatively the results of Ref.JHj]; the 
existence of weak and strong interacting TLSs in dis- 
ordered solids, and the corresponding strength of their 
interaction with the phonon field. As TLSs in KBr:CN 
were experimentally measured to have a coupling con- 
stant of 7 « 0.12 — 0.18eV with the phonon field, our 
calculations also verify that it is indeed the weak inter- 
acting t TLSs which are the relevant TLSs at low tem- 
peratures, dictating the universal behaviour. Thus, we 
are able to clearly identify the relevant TLSs in this par- 
ticular system. In Ref. [12j the plausibility that the same 
mechanism dictates universality in amorphous solids was 
argued for. The verification of this argument requires 
the identification of nearly inversion symmetric TLSs in 
amorphous solids, and the calculation of their coupling 
to the phonon field. 
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